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ABSTRACT : To investigate the genetic diversity of seven Cliinese indigenous meat goat breeds (Tibet goat, Guizhou white goat, 
Shannan white goat, Yichang white goat, Matou goat, Changjiangsanjiaozhou white goat and Anhui white goat), explain their genetic 
relationship and assess their integrity and degree of admixture, 302 individuals from these breeds and 42 Boer goats introduced from 
Africa as reference samples were genotyped for 11 microsatellite markers. Results indicated that the genetic diversity of Chinese 
indigenous meat goats was rich. The mean heterozygosity and the mean allelic richness (AR) for the 8 goat breeds varied from 0.697 to 
0.738 and 6.21 to 7.35, respectively. Structure analysis showed that Tibet goat breed was genetically distinct and was the first to separate 
and the other Chinese goats were then divided into two sub-clusters: Shannan white goat and Yichang white goat in one cluster; and 
Guizhou white goat, Matou goat, Changjiangsanjiaozhou white goat and Anhui white goat in the other cluster. This grouping pattern was 
further supported by clustering analysis and Principal component analysis. These results may provide a scientific basis for the 
characteristization, conservation and utilization of Chinese meat goats. (Key Words : Meat Goat, Genetic Diversity, Genetic 
Differentiation, Microsatellite) 



INTRODUCTION 

The domestic goat (Capra hircus) is one of the most 
important livestock species in animal husbandry. China 
raises the largest number of goats in the world and has a 
valuable resource of goat breeds. There are more than 20 
indigenous goat breeds in China, including goats for meat, 
wool, milk, mutton, and kid leather (Tu et al., 1989). Goats 
have a strong adaptability and are distributed in different 
environmental areas. There were several meat-type goat 
breeds in China. These meat goats, which are mainly 
distributed in the middle and southern regions of China, 
play an important role in the meat market. 

Genetic diversity and relationships among some 
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indigenous goat breeds in China have been studied using 
RFLP (Li et al., 1997; Zhao et al., 2011) and RAPD of 
genomic DNA (Li et al., 2000). In addition, there is some 
research on Chinese cashmere goats using SSR of genomic 
DNA (Di et al., 2010), but there is Umited information 
regarding Chinese indigenous meat goats based on 
microsatellite markers. 

In the past three decades, meat production performance 
has been the primary breeding objective. Genetic admixture 
among some breeds and inbreeding within the breeds is 
possible. This may have led to the loss of genetic structure 
and diversity in some of these breeds. Futhermore, some 
breeds have decreased rapidly in the number of sires and in 
population size such as seen in the shannan white and 
yichang white goats. Some breeds are even facing 
extinction. 

In this study, efforts were made to evaluate genetic 
diversity and relationships of 7 Chinese indigenous meat 
goat breeds and an African goat breed based on standard 
microsatellite markers. The goal was to have a better 
understanding of current diversity and the level of 
inbreeding in these breeds, and offer a contribution towards 
further improvement of conservation schemes. 
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MATERIALS AND METHODS 

Sample collection and DNA extraction 

Blood samples from 302 animals were collected from 
seven Chinese indigenous meat goat breeds (Table 1). 
Among these meat goat breeds, Tibet goats are dual- 
purpose meat and cashmere goats. The seven Chinese goat 
breeds come from seven different provinces, mainly 
distributed among the Yangtze River. Samples for the 
indigenous breeds were collected from their conservation 
farms. Additionally, we collected 42 Boer goats introduced 
from Africa as reference samples. Genomic DNA was 
extracted from all blood samples using a standard phenol- 
chloroform extraction method (Sambrook et al., 2001). 

Microsatellite genotyping 

Eleven microsatellites were used from the panel 
recommended by ISAG/FAO for goat diversity studies 
(Table 2). Forward primers were 5 '-labelled with 
fluorescent dyes (FAM or HEX). PGR amplifications were 
performed in 12 ^1 reaction volumes containing 50 ng of 
genomic DNA, 0.025 (jM of each primer, 250 (xM of each 
dNTP, 1.25 units of Taq polymerase and Ix Magnesium- 



Table 1. Sampling information of the 7 Chinese goat breeds 



Breed name and index Number Collection location 



Tibet goat (TB) 


45 


Changdu City, Tibet 


Guizhou white goat (GZ) 


38 


Zunyi City, Guizhou 


Shannan white goat (SN) 


48 


Ankang City, Shaanxi 


Yichang white goat (YC) 


44 


Yichang City, Hubei 


Matou goat (MT) 


42 


Changde City, Hunan 


Changjiangsanjiaozhou 


46 


Nantong City, Jiangsu 


white goat (CSJ) 






Anhui white goat (AH) 


39 


Feidong county, Anhui 



free PGR buffer (Takara, Japan). Amplifications were 
carried out using the GeneAmp PGR 9700 thermocycler 
(AppUed Biosystems), with the following cycUng 
parameters: 94°G for 5 min followed by 32 cycles of 94°G 
for 30 s, annealing at 53 to 60°G for 30 s, 72°G for 30 s, and 
a final step at 72°G for 8 min. PGR products were diluted 
1/3 to 1/4 and 0.75 |li1 of each diluted product was then 
mixed with an internal standard (GeneScan™, 500 LIZ ™, 
Applied Biosystems) according to the manufacturer's 
instructions. Genotyping was carried out on an ABI 3130x1 
automated capillary sequencer. 



Table 2. Information of 11 microsatellite loci used in this study 



Locus 


Primer sequence 

(5'-3') 


Forward 
label 


Set 


Tm 

(°C) 


At 


Rt 


McM527 


GTCCATTGCCTCAAATCAATTC 
AAACCACTTGACTACTCCCCAA 


HEX 


1 


55 


8 


6.5 


ILSTS005 


GGAAGCAATTGAAATCTATAGCC 

TGTTCTGTGAGTTTGTAAGC 


FAM 


1 


60 


7 


6.5 


SRCRSP9 


AGAGGATCTGGAAATGGAATC 
GCACTCTTTTCAGCCCTAATG 


FAM 


1 


58 


7 


4.5 


OarFCB20 


GGAAAACCCCCATATATACCTATAC 

AAATGTGTTTAAGATTCCATACATGTG 


FAM 


2 


58 


9 


7.9 


SRCRSP5 


GGACTCTACCAACTGAGCTACAAG 
TGAAATGAAGCTAAAGCAATGC 


FAM 


2 


53 


13 


8.0 


TGLA53 


GCTTTCAGAAATAGTTTGCATTCA 
ATCTTCACATGATATTACAGCAGA 


HEX 


2 


55 


11 


7.6 


ILSTSOU 


GCTTGCTACATGGAAAGTGC 
CTAAAATGCAGAGCCCTCCC 


HEX 


2 


53 


11 


6.6 


SRCRSP23 


TGAACGGGTAAAGATGTG 
TGTTTTTAATGGCTGAGTAG 


FAM 


3 


58 


9 


8.1 


SPS1I3 


CCTCCACACAGGCTTCTCTGACTT 
CCTAACTTGCTTGAGTTATTGCCC 


FAM 


3 


60 


10 


8.1 


SRCRSP7 


TCTCAGCACCTTAATTGCTCT 
GGTCAACACTCCAATGGTGAG 


HEX 


3 


55 


8 


6.0 


INRA063 


GACCACAAAGGGATTTGCACAAGC 
AAACCACAGAAATGCTTGGAAG 


HEX 


3 


58 


13 


8.5 



Tm = Annealing temperature; At = Total number of alleles per locus; Rt = Allelic richness over all samples. 

Hoffmann, I., P. A. Marsan, S. F. Barker, E. G Cothran, O. Hanotte, J. A. Lenstra, D. Milan, S. Weigend and H. Simianer. 2004. New MoDAD marker sets 
to be used in diversity studies for the major farm animal species: recommendations of a joint ISAG/FAO working group. In: 29th International 
Conference on Animal Genetics. FAO, Meiji University, Tokyo, Japan. 
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Statistical analysis 

The absence of linkage between loci was verified using 
GENEPOP 3.4 (Raymond and Rousset, 1995). The 
polymorphism parameters were obtained using Excel 
Microsatellite TooUdt Version 3.1 (http://animalgenomics. 
ucd.ie/sdepark/ms-toolkit/index.php) and POPGENE 1.31 
(Yeh and Boyle, 1997). The F-statistic values (FIT, FIS, and 
FST; Weir and Cockerham, 1984), together with the allelic 
richness (AR) of each breed, were estimated with FSTAT 
2.9.3 (Goudet, 2001). Nei's standard genetic distance Ds 
(Nei, 1972) and Nei's genetic distance (Nei et al., 1983) 
between breeds were calculated using the DISPAN package. 
Neighbour-joining (NJ) dendrograms and the Unweighted 
pair group method with arithmetic mean (UPGMA) were 
constructed based on and Dg. The robustness of the 
dendrograms was evaluated using a Bootstrap test of 1,000 
resamplings of loci with replacement. Principal component 
analysis (PCA) was performed to reveal major patterns of 
genetic variabiUty and clustering of breeds based on allele 
frequencies using MVSP 3.1 (http://www.kovcomp.com). 
Based on the genotypes at the 11 microsatellite markers, 
individual animals were clustered into a given number of 
populations and assigned probabilistically to clusters 
inferred with a Bayesian approach implemented by the 
STRUCTURE software (Pritchard et al., 2000). 

RESULTS AND DISCUSSION 

Polymorphism of microsatellite loci 

Results showed that 11 loci were polymorphic. A total 
of 106 alleles were detected, and the mean number of 
alleles across the 11 loci was 9.6. The number of alleles per 
locus (At) ranged from 7 for ILSTS005 and SRCRSP9 to 13 



for 1NRA063 and SRCRSP5, which is greater than that 
reported by FAO-ISAG (Table 2). Allelic richness over all 
samples per locus (Rt) was measured at between 4.5 
(SRCRSP9) and 8.5 (INRA063) and the mean allelic 
richness across 11 loci in our samples was 7.1. Generally, 
these loci were polymorphic, and the polymorphism 
information content (PIC) across the 11 loci ranged 
between 0.545 (SRCRSP7) and 0.778 (1NRA063). 

The HWE was tested for all breed-locus combinations. 
Significant (p<0.05) deviations from a HWE were observed 
for 13 (14.8%) of 88 breed-locus combinations (data not 
shown). But no significant (p>0.05) deviation from a HWE 
was detected for either a single locus across all breeds or a 
single breed across all loci. 

Genetic diversity 

There was a high allelic variation found in the Chinese 
goat breeds for 11 loci. A summary of the identified 
polymorphisms from all of the meat goat breeds we tested 
is hsted in Table 3. For the Chinese goats, the observed 
heterozygosity (Ho) ranged from 0.697 to 0.738 with the 
lowest value found in Tibet goats and the highest in Matou 
goats. The outgroup Boer goat breed was 0.634 (SD = 
0.024). 

There were notable differences in the H and AR among 
these breeds (Table 3). The unbiased gene diversity (H) 
values varied between 0.738 in Tibet goats and 0.781 in 
Matou goats (mean value: 0.762; SD: 0.014). Allelic 
richness (AR) was the highest in the Anhui breed (AR = 
7.35) and lowest in Tibet (AR = 6.21), with a mean of 6.78. 
Allelic richness (AR) in the Boer was 6.07. For Chinese 
goat breeds, the mean number of alleles (MNA) among 
these goat breeds ranges from 6.27 (Tibet) to 7.36 (Anhui). 



Table 3. Sampling information and basic parameters for the genetic diversity associated with the 8 meat goat breeds analyzed in this 
study 



Breeds and codes 




Allelic diversity 




Genetic diversity 


FIS 


TNA 


MNA 


NEA 


NPA 


AR 


H(SD) 


Ho(SD) 


Tibet goat (TB) 


69 


621 


3.93 


2 


6.21 


0.738(0.013) 


0.697(0.022) 


-0.056* 


Guizhou white goat (GZ) 


75 


6.82 


4.25 


3 


6.82 


0.764(0.010) 


0.711(0.020) 


0.071** 


Shannan white goat (SN) 


11 


7.00 


4.32 


4 


6.88 


0.759(0.016) 


0.706(0.021) 


0.069** 


Yichang white goat (YC) 


73 


6.64 


4.20 


8 


6.59 


0.748(0.011) 


0.711(0.016) 


-0.051** 


Matou goat (MT) 


11 


7.00 


4.55 


4 


6.94 


0.781(0.014) 


0.738(0.020) 


0.036 


Changjiangsanjiaozhou 


lA 


6.73 


4.49 


2 


6.64 


0.774(0.014) 


0.723(0.016) 


0.066*** 


white goat (CSJ) 


















Anhui white goat (AH) 


81 


7.36 


4.59 


5 


7.35 


0.771(0.011) 


0.718(0.021) 


0.063 


Mean 










6.78 


0.762(0.014) 


0.715(0.017) 




Boer goat (BR) 


67 


6.09 


4.15 


3 


6.07 


0.734(0.013) 


0.634(0.024) 


-0.137*** 


ALL Mean 










6.69 


0.759(0.013) 


0.705(0.018) 





TNA = Total imrnber of alleles; NEA = Number of effective alleles; MNA = Mean number of alleles; AR = Allelic richness; NPA = Number of private 
alleles; H = Unbiased gene diversity; Ho = Observed heterozygosity; SD =Standard deviation; HWE = Number of loci that deviate from a Hardy 
Weinberg equilibrium; FIS = Population inbreeding coefficient. *** p<0.001; ** p<0.01; * p<0.05. 
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A total of 31 private alleles (NPA) were amplified in our 
cohort, 28 unique to Chinese goats and three unique to Boer 
goats. The NPA of the Yichang goats was particularly high 
(NPA = 8), representing 29% of the total NPA. Most of the 
private alleles (26) were at very low frequencies of below 
5%. Four unique to Chinese goats and one to a Boer goat 
showed a frequency that exceeded 5%. 

The inbreeding coefficients (FIS) for five Chinese 
breeds were positive, and three of them (GZ, SN and CSJ) 
reached significant levels (p<0.01. Table 3). The global FST 
value of 0.056 for Chinese meat goats in this study was 
similar to that of the Chinese goats (0.063, Di et al., 2011; 
0.054, Li and Valentini, 2004) and the European and Middle 
Eastern goats (0.069, Canon et al., 2006), but lower than the 
value of 0.143 of the Southeast Asian goats (Barker et al., 
2001). 

The relative small number of microsatellites (11 loci) 
may bring some limitations to the analysis. But in this 
article, the polymorphism of selected microsatellites was 
rich and these goat breeds had a relative large sample size. 
These could reduce the limitations caused by the relative 
small locus number. 

Genetic phylogeny 

Genetic phylogeny was analysed with three methods 
following that of Peter et al. (2007). Using STRUCTURE 
2.1, we determined an optimum K value at 4 following 
Evanno et al. (2005), who found that an ad hoc quantity 
based on the second order rate of change of the likelihood 
function with respect to K (AK) did show a clear peak at the 
true value of K. Figure 1 shows a clear clustering (K = 4), 
which was in accordance with breed types. Tibet goat, 

<b A .O ^ -i" ^ 
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Figure 1. Population structure of eight meat goat breeds displayed 
with individual (top half of figure) and population (bottom half of 
figure) Q-matrices at K = 4 using STRUCTURE 2.1 following ten 
independent runs for 2<K<6 with 10* Markov Chain Monte Carlo 
repetitions. TB = Tibet goat; GZ = Guizhou white goat; SN = 
Shannan white goat; YC = Yichang white goat; MT = Matou goat; 
CSJ = Changjiangsanjiaozhou white goat; AH = Anhui white goat; 
BR = Boer goat. 



56 



50 


81 



















CSJ 
AH 
MT 
GZ 
YC 
SN 
TB 
BR 



Figure 2. Unweighted pair group method with arithmetic mean 
(UPGMA) tree based on Nei's DA distances for 8 meat goat 
breeds. TB = Tibet goat; GZ = Guizhou white goat; SN = Shannan 
white goat; YC = Yichang white goat; MT = Matou goat; CSJ = 
Changjiangsanjiaozhou white goat; AH = Anhui white goat; BR = 
Boer goat. 

which belongs to the dual-purpose meat and cashmere goats, 
was genetically distinct and was the first to separate. The 
other Chinese goats were then divided into two sub- 
clusters: SN and YC in one cluster; and GZ, MT, CSJ and 
AH in the other cluster 

In addition, an UPGMA tree was constructed on the 
basis of the genetic distances with relatively high 
bootstrap values (Figure 2). As expected, the BR was most 
distinct for other breeds first, and the TB was divided from 
other Chinese goat breeds. There were two clusters found 
among the remaining Chinese goats. This grouping pattern 
was further supported by PCA analysis. Since phylogenetic 
reconstruction may not easily take into account the effects 
of admixture between breeds, we performed PCA method to 
further investigate possible genetic relationships between 
Chinese goat breeds. The first principal component (PC) 
explains 31.2% of the observed genetic variation, and the 
second and third PCs resolve 23.2% and 16.0% of this 
variation, respectively. A two-dimensional scatter plot for 
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Figure 3. Principal component analysis of 7 Chinese indigenous 
meat goat breeds. Axis 1 and Axis 2 represent the first two 
principal factors, respectively. TB = Tibet goat; GZ = Guizhou 
white goat; SN = Shannan white goat; YC = Yichang white goat; 
MT = Matou goat; CSJ = Changjiangsanjiaozhou white goat; 
AH = Anhui white goat; BR = Boer goat. 
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the Chinese goat breeds we analysed (Figure 3) reveals that 
they are clustered into several different groups and that 
reflects the geographical distribution of these breeds, with 
the secend component corresponding to distributions from 
east to west, respectively. 

Of the Chinese indigenous goats, the Tibetan goat was 
the most distinct according to structure plot (Figure 1). An 
early study proposed the Tibetan Plateau to be a place of 
origin for Chinese goats (Chen et al., 2006). Combined with 
earHer studies (Chen et al., 2005; Naderi et al., 2008; Wu et 
al., 2009) and archaeological records, we hypothesize that 
the Tibetan goat was an ancient breed and might carry 
unique adaptation alleles. Therefore this breed should be 
considered as a priority group for conservation and further 
research. 

According to our investigation, the individual 
exchanges among meat goat breeds for improving meat 
quantity (Pan et al., 2009; Mao et al., 2010) might have led 
to a close relationship among meat goat breeds (Figure 1). 
To maintain the present genetic diversity and structure of 
these breeds, genetic exchanges between breeds must be 
controlled on every conservation farm. 
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